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I suggest simple method for the search of Dark Matter particles and some related particles which 
allows to measure reliably their mass and spin in a wide class of models for Dark Matter. 
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I. INTRODUCTION 

About 25% of the Universe is made from Dark Mat- 
ter (DM). There are a number of models in which DM 
consists from particles similar to those in SM with addi- 
tional discrete quantum number, which I denote here as 
D-parity. For known particles D = 1, for new DM parti- 
cles (DMP) D = — 1, and D-parity conservation ensures 
stability of lightest from these new particles. We will con- 
sider such models for DM, in which there is more than 
one particle with parity D = —1, in particular, neutral 
DMP, denoted below as D, and its more heavy lightest 
charged counterpart Z? ± . It is assumed that all these 
D-particles have identical spin sd (1/2 or 0). 

• The well known examples of such models for DM are 
given by MSSM and NMSSM. Here D-parity is another 
name for R-parity, D is neutralino and is chargino, 
here Sd = 1/2. There exists a vast literature here, see 
for example, [l|. 

• The second example is given by inert doublet model 
(IDM) @. In notations 1 31 , this model contains one stan- 
dard Higgs doublet cj>s, responsible for electroweak sym- 
metry breaking and generation of fermions and gauge 
bosons masses as in the Standard Model (SM), and an- 
other scalar doublet (fir), which doesn't receive vacuum 
expectation value and doesn't couple to fermions. Four 
degrees of freedom of the Higgs doublet (fi$ are as in the 
SM: three Goldstone modes and one mode becomes the 
Higgs boson h. All the components of the second scalar 
doublet (ftp are realized as massive scalar D-particles: 
two charged and two neutral ones D&nAD A . By con- 
struction, they possess a conserved multiplicative quan- 
tum number - D-parity, the lightest particle among them 
is considered as a candidate for DM particle. In our no- 
tations the lightest D-particle is D. Here sd = 0. 

Detailed limitations for masses of DM-particles from 
cosmology and collider physics are discussed in many 
papers (see e.g. [l|, Q). Below we have in mind 
M D < 70 GcV and M D ± > 80 GeV. 

We consider production of such D-particles at e + e~ 
linear collider (e.g. at ILC) with electron beam energy 
E e = y/s/2, provided that the family of D-particles with 
smaller weight than E e is settled by D, D ± and no 
more than one neutral (for example D A ). All processes 
below we treat as basic reactions, e + e~ — > D + D~ or 
e+e~ — > DD A with subsequent decay of D ± or D A . 

The problems which should be solved in the discussed 



experiments are following. 

(1) To observe process with production of D-particles. 

(2) To evaluate mass of DMP Md and masses of other 
D-particles. 

(3) To evaluate spin of D-particles Sd- 

(4) To obtain some additional information about inter- 
actions of D-particles. 



II. TWO TYPE OF D-PARTICLES, D AND D . 

In this section we consider the case when the set of D- 
particles in the energy region of ILC is settled by and 
D. In the considered models the only interaction that 
couples D and D ± with SM matter is gauge interaction 



D ± DW T . 



(1) 



Therefore, D ± — > DW^ with either on shell or off shell 
VF's is the only decay channel for D^. 

A. Production, decay, signature 

We suggest to use D + D _ production at ILC: 

e+e- -> D+D- -> DDW + W~ . (2) 

The lab system here id c.m.s. for e + e~. In this system 
energies, 7-factors and velocities of D ± are 



E± =E e , j±=E e /M D ±, P± 



\-MlJEl (3) 



Neglecting quantity 1/4 — sin 2 8^ (describing 7 — Z 
interference) the cross section of this process reads 
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numerical factor). These cross sections are of the same 
order as a(e + e~ — > fj.fi) (which is given by dU at (3 = 1). 
That are huge cross sections for ILC. 

• If Mp± — Mn > Mw ■ the produced Z) ± decays to 
DW ± with mass shell W only. We suggest to observe 
the following final state systems: 

• Two dijets from qq decay of W + and W~ , with ef- 
fective mass Mw- For this channel the cross section is 
- [0.676 2 « 0.45] • a(e+e- -> D+D~). 

• One dijet from qq decay of W + or W~ plus (i or e 
from fiv or ev or tv — > fivvv or tv — > evvv decay 
of W~ or W + . For this channel the cross section is 
[2 • 0.676 • 2 • (1 + 0.17)0.108 » 0.33] • o{e+er -> D+D~) 
(here 0.17 is a fraction of /j or e from the decay of r). 

Typical event will have large missing energy carried out 
by neutral and stable Z?-particles, the observed momenta 
are strongly a-collinear, the missing momentum is large 
and its direction is far from e + e~ collision axis. 

The background final state with the same kinematics 
can be produced in SM only if missed energy and p± 
are carried off by neutrino(s). The value of correspond- 
ing cross section is at least by one electroweak coupling 
squared g 2 /47r or g' 2 /4tt less than <j(e + e _ — > fifi) with 
g 2 /4:Tr ~ g' 2 /47r ~ a. Therefore, cross sections for such 
background processes are about two orders less than the 
cross section of the process under discussion. The same 
estimate is valid for all reactions, considered below. 

• If Mni — Mo < M\v, the only mode for charged 
decays ifQ DW* + . The structure of the final state and 
the value of cross sections are the same as in previous case 
with a single difference - dijet effective mass M* is now 
not peaked around resonance value Mw but distributed 
in some interval below M D ± — Mo- The form of this 
distribution depends on the spin of D-particles So- 

• The signature for the process in both cases is: 



energy and momentum of W in this frame are 



Two dijets or one dijet plus e or /i with large 
missing energy and large a-collinearity, with 
cross section ~ er(e+e~ — > fi + fi~). Typically 
these dijets (or dijet and lepton) move in the 
opposite hemispheres. 



(6) 



B. Parameters of and D 



We denote 



A(s; si, s 2 ) = Y s 2 + sj + ,s 2 .-2ssi-2ss 2 -2siS2- (7) 

The case Mo± > Moo + Mw ■ In the rest frame of 
D ± we deal with 2- particle decay D ± -> DW ± . The 



1 In the description of reaction we denote by W and Z both W and 
W* , Z and Z* . Here W* and Z* mean qq dijet or dilepton system 
with quantum numbers of W or Z but with lower effective mass. 
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The effective mass of dijet is My/ with accuracy to W 
width. 

Denoting W escape angle in D + rest frame relative to 
the direction of D + motion in the lab system by 9 and 
c = cos 9 we have energy of W + in the lab system 

E w =-f±(E w + cf3 ±P r ). 

Therefore, energies of dijets W's are distributed within 
the interval 

(£(-)=7±(^-/3±P r ), E (+) = 1± {E w +p ± f )). (9) 

Masses. The end point values -E(±) give two equations 
for evaluation of masses D ± andD . 
In particular, it is useful to note that 

=7^4' + b'') 2 , E ( _ )+ E (+) = 2 1± E W . (10) 
Therefore at large enough electron energy (at 7 ^> 1), 
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At finite 7 the exact equations are more complex. 

The accuracy of this procedure is determined by the 
accuracy of measurement of dijet energy together with 
its effective mass and by a width of Z? ± (if the latter is 
large). In particular, at so = the decay — > DW^ 
width is 



r\3 
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The r/M D ± ratio is below 0.1 at M D ± < 500 GeV. 

The distribution of these dijets in energy is uniform, 
dN(E) oc dE since there is no correlation between escape 
angle of W in the rest frame of £) ± and production angle 
of Z) ± . When the width of £> ± is not small, this dis- 
tribution become non-uniform near the end points. The 
measuring of fine structure of this distribution near the 
end point will give the total width at least roughly. 

• If Mo± > Moo + Mw , the only decay channel is 
D + —> DW* + . All discussed above results are valid for 
each separate value M*, with evident change in all equa- 
tions M w -t M*. 

The energy and M* distributions for each pair of dijets 
are independent from each other. 

The masses of M D ± and Mo are evaluated in this case 
for each measured value of M* even with the best accu- 
racy than in the previous case since in this case proper 
width of Z) ± is low enough. 

Spin. After evaluation of M D ±, the cross section of 
e + e~ — > D + D~ process is calculated precisely in QED 
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for each sd- The cross section for sd = ([5]) is more 
than four times less than that for Sjj = 1/2 (JD). This 
big difference allows conclude definitely about spin of D- 
particles. 

Other properties for sd = 1/2. If D ± are fermions, 
their spins are correlated to each other and to longitu- 
dinal polarization of collided electrons. Observation of 
process e+e" -> D+D- -> D°D°jj£ + u'b allows to find 
the sign of charge of dijet W = qq in each separate case. 
It allows to study the charge and polarization asymme- 
tries for accessing of more detail properties of D-particles 
(e.g. ratio of D + D~j to D + D~Z couplings). 



III. TWO NEUTRAL D-PARTICLES 

Here we consider the case when the set of D-particles 
in the energy region of ILC consists of single charged D- 
particle D and two neutral D-particles with identical 
CP, D and £>i with M D ± > M Dl > M D , where D can 
be DMP while D\ can disappear during lifetime of the 
Universe. In the discussion below we neglect an exotic 
case of approximate degeneracy between D and D\. 

• Let the only interaction that couples D or D\ and 
Z? ± with SM matter is the standard gauge interaction 
(|T|) . In this case process ^ is supplemented by similar 
processes with one or two D± instead of D in the final 
state: 



D + D~ 



DDW+W- 
DD X W+W- 
D X D X W + W' 



(a), 
(b) 
(c) 



(13) 



with the same signature ^ . The probability of the decay 
Z) ± — > £>iW ± is lower than the probability of the decay 
Z) ± — > DW due to smaller final phase space volume in 
the first case. Therefore a(a) > a(b) > cr(c). 

The description of basic process e + e~ — > D + D~ 
coincides with that presented above. Each decay 
D ± -> DW ± and D ± -> D 1 W ± is described as 
above, the only new point is that the end points Et-) 
and -EV+) for energy distribution of W's from decay 
Z? ± — > DW ± are given by equations ([9]) while the end 
points E}s and Eh-, for energy distribution of W's from 

decay D ± — » £>iW ± are given by the same equations 
but with the change Mr> to Mu x - It is evident that 
J5(_) < E}-. < Eh-, < E( + y Therefore in this case the 
same procedure as above allows to obtain masses M D ± 
and M D . 

To evaluate Mr> 1 let us remind that for each type of 
decay the energy distribution of W in the lab. system 
is uniform. The energy distribution of W in the lab. 
system is the sum of two uniform distributions with the 
described above end points. The end points Ely Eh-. 
are marked by the steps in the density of event energy 
distribution. They can be used for evaluation of Mu ± . 

The distribution of dijets in the effective mass can be 
different. If all masses are peaked around My/, it means 



■D+D- 



that M D ± -M D > M D ± - M Dl > M w . If there are 
dijets with effective mass Mw and those with lower ef- 
fective mass, the former appear from Z) ± — > DW ± decay 
and the latter - from D ± — > DiW^ decay. They can 
be used for evaluation of Mu and Md 1 separately. If 
the effective masses of all dijets are below Mw, we have 
Mw > M D ± — Md > M D ± — Mdj. In this case men- 
tioned steps in the dijet energy distribution at each M* 
will be added by steps in the distribution in M* for eval- 
uation of Mo and Mr> 1 . 

Both channels of decay must be taken into account 
at comparison of the measured cross sections with (HJ) , ([5]) 
for evaluation of spin sp. 

• In some models interaction D\Dh is allowed ad- 
ditionally. We discuss this opportunity assuming that h 
is anticipated Higgs boson of SM with mass about 120 
GeV, very small width and dominant bb decay channel. 
Modifications of description for another properties of h, 
observed at LHC before operations at LC, is simple. In 
this case reactions (TTSTb)) and ([TBT c')) are modified as 

DLW+W- -^DDW+W-h, (b) 
DxDxW+W- -^DDW+W-hh. (c) 

(14) 

The signature @ is naturally added by one or two bb 
dijets from the Higgs decay. Using of b tagging allows to 
distinguish these dijets from W dijets which is necessary 
for the evaluation of masses. If the ratio (Md x —Mo)/Mh 
is small enough, the D\ time of life can be so large that 
one can observe it in the detector. 



IV. TWO NEUTRAL Z)- PARTICLES WITH 
OPPOSITE CP-PARITY 

In the IDM together with D and D ± the one more neu- 
tral scalar particle D A exists (with mass M d a > Md) 
3]. Here CP parities of D or D A cannot be defined sep- 
arately since they do not interact to fermions, but their 
relative parity is fixed, they have opposite CP-parities. 
It results in the gauge interaction D A DZ and instabil- 
ity D A via the decay D A —> DZ (in addition to gauge 
interaction D A D + W~ similar to that for D, DD + W~). 
Similar particles, interaction and production channel are 
possible also in some models for fermionic Dark Matter. 

So, we discuss now the case when additional neutral 
D-particle D A . In this discussion we neglect case of 
approximate degeneracy between D and D A . 

• If M d a + M d < 2E e , M d a < M D ±, The most im- 
portant new process is 

e+e- -> DD A -> DDZ. (15) 

Instead of ©, D A energy, 7- factor and velocity are 
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In the IDM (at so = 0) the cross section 

E A 



a(e + e~ -> DD ) = 



3(s-M|) 2 ^ 2E e -E A 



(17) 



This cross section is of the same order of value as 

The signature of this process is similar to that given 
by eq. ([H]) : 



One qq dijet or e + e~ or [V pT pair with identi- 
cal effective mass and energy distributions and 
with large missing E±. 



(18) 



If M d a — Md > Mz, observable final state is Z, which 
is seen as hadronic dijet or e + e~ or with effec- 

tive mass equal to Mz- Energy distribution of this Z is 
given by equations similar to (O,©. End points of this 
distribution allow to evaluate masses Mo and Mo a . 

If M d a —Mo < Mz, observable final state is Z* , which 
is seen as hadronic dijet or e + e~ or with identical 

spectra of effective mass. For each value of this effective 
mass energy distribution of this Z* is given by equations 
similar to End points of this distribution allow 

to evaluate masses Mo and M d a . 

After the study of process (|15p one must to study pro- 
cess ([5]) and cascade reactions 



>D + D~ 



DDW+W- , 
D A W ± DW T - 
D A W+D A W~ 



(a) 



■DDW+W-Z, (b) 
-^DDW+W-ZZ. (c) 
(19) 

The decay D+ — > D A W ± is described by the same equa- 
tion as the decay D+ — > DW+. Its probability is lower 
than that for decay D+ —> DW+ due to smaller final 
phase space volume. Therefore, a (a) > a{b) > cr(c). 

The signature of the process (fl9lf a)). just as the process 
(TT97 b). (c)) for invisible decays of Z, is given by ©■ Just 
as in the analysis of sect. IIII1 the study of end points in 
the energy distribution of W allows to obtain Mo and 
Mo± and steps in this distribution can be used to verify 
value of M d a , obtained via a study of process with the 
signature (TTg)) . 

The signature of processes (fPW b)) and (fl9T c)) for vis- 
ible decays of Z is similar to §§§ with adding of dijet or 
£ + £~ pairs which represent Z or Z*. If one can distin- 
guish dijets from Z and those from W, the energy distri- 
bution of W in these processes can be used to enhance 
data massive for evaluation of masses, discussed above. 

Summation of these 3 cross sections gives total cross 
section of e + e~ — > D + D~ process, which value gives us 
the value of spin s o ■ 

• If M d a + Md < 2E e , M d a > M d l , the analysis of 
sect. [IT] is valid completely for the final states with signa- 
ture (O - reaction e + e~ — > D + D~ . 



The second series of processes is e+er -> DD A with 
two channels of D decay and different signatures 



DD 



DDZ , 
DD ± W T -> DDW+W- 



(a) 
(b) 



(20) 



The process pOTa)) is the process (TT5|) . It can be 
analysed just as it was discussed earlier. The process 
([2"0Tb)) is a cascade process. It can be eliminated 
from mass analysis of process ([2]) by using the fact 
that in difference with the process the observable 
decay products of this process move typically in one 
hemisphere. 



V. SUMMARY 

We present simple and robust method for discovery 
candidates for DM particles and evaluation of its mass 
and spin (with evaluation of mass and spin of some other 
particles with the same -D-parity) at ILC/CLIC. 

These particles will be discovered via observation of 
processes with signature ([6]) , ([TBI) and with cross section 



of the order of a(e + e 



which is huge for LC. 



The masses of these particles will be obtained via mea- 
suring the details of energy distribution of dijets (repre- 
senting W+ or Z) near the end points. The cross section 
measurement of process with signature ^ , (fTSJ) and sim- 
ilar signature for the derivative processes with cascade 
decay allows to determine the spin so of considered can- 
didate for DM particle by comparison with simple SM 
calculation (the cross sections for so = 1/2 is approxi- 
mately 4 times larger than that for = 0) . 

One of proposed processes, e + e~ — > D + D~ , was con- 
sidered earlier in respect of discovery of neutralino as 
DMP, etc. (see e.g. @). However, I never saw such ap- 
proach for simultaneous evaluation of masses and spins 
of D-particles irrespective to details of model. The ad- 
vantages of presented approach are following: 

1. The cross section of each suggested process is substan- 
tial part of the total cross section ofe + e~ annihilation at 
considered energy (typically up to tens percents). 

2. The signature is clear. 

3. The background is very small ( typically, ~ 1 % from 
the observable effect). 

4- The kinematics is very simple, allowing reliable ex- 
traction of quantities under interest from the data. 
We have no idea how one use similar approach at 
hadronic colliders. 

I am thankful D.Yu. Ivanov and K.A. Kanishev for 
discussions. The work was supported by grants RFBR 
08-02-00334-a, NSh-3810.2010.2 and Program of Dept. of 
Phys. Sc. RAS "Experimental and theoretical studies of 
fundamental interactions related to LHC." 
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I. INTRODUCTION 

About 25% of the Universe is made from Dark Mat- 
ter (DM). There are a number of models in which DM 
consists of particles similar to those in SM with addi- 
tional discrete quantum number, which I denote here as 
D-parity. For known particles D = 1, for DM particles 
(DMP) D = — 1, and D-parity conservation ensures sta- 
bility of the lightest particle with D = — 1. We consider 
such models for DM, in which there is more than one D- 
odd particle, in particular, neutral DMP D and its more 
heavy lightest charged counterpart D ± and all these D- 
particles have identical spin (1/2 or 0). 

• The well known examples of such models for DM are 
given by MSSM and NMSSM. Here D-parity is another 
name for R-parity, D is neutralino and D is chargino, 
here sd = 1/2 (see for example [l|). 

• The second example is given by inert doublet model 
(IDM) [2]. In notations d, this model contains one stan- 
dard Higgs doublet ips, responsible for electroweak sym- 
metry breaking and generation of fermions and gauge 
bosons masses as in the Standard Model (SM), and an- 
other scalar doublet <f>jj, which doesn't receive vacuum 
expectation value and doesn't couple to fermions. Four 
degrees of freedom of the Higgs doublet cj>s are as in the 
SM: three Goldstone modes and one mode which be- 
comes the Higgs boson h. All the components of the 
second scalar doublet <pjj are realized as massive scalar 
D-particles: two charged D ± and two neutral ones D 
and D A . By construction, they possess a conserved mul- 
tiplicative quantum number - D-parity, the lightest par- 
ticle among them is considered as a candidate for DM 
particle, here Sd = 0. 

Limitations for masses of DM-particles from cosmology 
and collider physics are discussed in many papers (see e.g. 
[l|, 0| )• These limitations allow existence of discussed 
particles with masses lower than electron beam energy 
of ILC/CLIC. We don't discuss here the case of mass 
difference M D ± — Mr> or M D ± — M d a lower than about 
10 GeV. 

In the considered models the interaction of D-particles 
with SM particles appears only from the covariant deriva- 
tive in the kinetic term of the Lagrangian, that are gauge 
interactions 



with standard electroweak gauge couplings g and g'. 

We consider production of such D-particles at ILC 
with electron beam energy E e = \/~s/2, provided that 
the family of D-particles with smaller mass than E e is 
confined to D, D^ and no more than one additional neu- 
tral D A . All processes below we treat as basic reactions 



(<*) 



e ' e 



D+D- 



(b) 



(2) 



with subsequent decay of D^- or D A . 

The problems which should be solved in these experi- 
ments are the following. 

(1) To observe unambiguously the process with the pro- 
duction of D-particles. 

(2) To evaluate the mass of DMP Mn and masses of other 
D-particles. 

(3) To evaluate the spin of D-particles sd- 

(4) To obtain some additional information about inter- 
actions of D-particles. 

In all presented cross sections we neglect quantity 1 / 4— 
sin 2 9w, describing j — Z interference. 



II. TWO TYPES OF D-PARTICLES, D AND D ± . 

In this section we consider the case when the set of 
D-particles in the energy range of ILC is confined to D ± 
and D with only decay channel D ± — > DW ± (with either 
on shell or off shell W's). 



A. Production, decay, signature 

We suggest to use reaction ((2^) with decay 
D ± -> dW ± : 



D+D- 



DDW + W~ 



(3) 



D+D-7, D+D-Z, 
D+DW~, D+D A W~, D A DZ 



(1) 



1 The interactions ((TJ form complete set of non-diagonal interactions 
of .D-particles. Therefore, the probability of discussed decay is 1 
independent on what is observed in the final state W or W* , where 
W* is of shell state of W, i.e. qq jets or Iv system with quantum 
numbers of W and lower effective mass M* . In the lists of reactions 
below we don't distinguish W and W* Z and Z* . 
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In the lab system (coincident with c.m.s. for e + e ) en- 
ergies, 7- factors and velocities of D ± are 



E± = E e , ~/±=E e /M D ±, p± 



l-M 2 /E 2 . (4) 



The cross section of this process reads 



cr(e+e- -> D + D~) = 
— /3(3-/? 2 )(l+ ,/ 1|<f2 , 2 ) at s D = -; (5) 



(s-M 2 ) 2 



33 I 1 i % 



3s 



■/3 d 1 + 



( S -M 2 ) 2 



at sn = 0; 



(6) 



R 



(1/2) 



(o) _ cot 2 (20 w ) 



16sin 4 (2<V) 



K z — 



4sin 2 (26»vy) 



(In different models for DMP the relative value of Z con- 
tributions can differ from these values by a simple 
numerical factor). These cross sections are of the same 
order as a(e + e~ — > fxfj,) (which is given by (© at /3 = 1). 
For ILC these cross sections are huge. 

If Md± — Mo > Mw the produced D ± decays to 
DW^ with on-mass shell W only. We suggest to ob- 
serve the following final state systems: 

• Two dijets from qq decay of W + and W~, with ef- 
fective mass Mw- For this channel the cross section is 
[0.676 2 « 0.45] • a(e + e" -> D+D~). 

• One dijet from qq decay of W + or W~ plus \i or e 
from [iv or ev or tv —¥ fivvv or tv — > evvv decay 
of W~ or W + . For this channel the cross section is 
[2 • 0.676 • 2 ■ (1 + 0.17)0.108 w 0.33] • cr(e + e- £>+£>-) 
(here 0.17 is a fraction of or e from the decay of r). 

Typical event will have the large missing transverse 
energy Ej_ carried out by neutral and stable particles 
D. The background is given by SM processes with the 
same kinematics and with large missed transverse energy 
E±, carried off by neutrino(s). The value of correspond- 
ing cross section is at least by one electroweak coupling 
squared g 2 /An or g' 2 /4tt less than <j(e + e~ — > /i/i) with 
g 2 /4ir ~ g' 2 /47r ~ a. Therefore, cross sections for such 
background processes are by about two orders less than 
the cross section of the process under discussion. The 
same estimate is valid for all reactions considered below. 

If M D i — Mq < Mw, the only mode for charged D^ 1 
decays is DW* + . The structure of the final state and the 
value of cross sections are the same as in the previous case 
with only difference - dijet effective mass M* is now not 
peaked around resonance value Mw but it is distributed 
in some interval below M D ± — Mo- The form of this 
distribution depends on the spin of D-particles so- 

• The signature for the process in both cases is: 



Two dijets one dijet plus e or /x with large 
missing energy and large a-collinearity + noth- 
ing, with cross section ~ cr(e + e _ — > /i+^ _ ). 
Typically these dijets (or dijet and lepton) 
move in the opposite hemispheres. 



(7) 



B. Parameters of and D 



We denote 



A(s;si,s 2 ) = ys 2 + s 2 + s 2 -2ssi-2ss 2 -2siS2. (8) 

The case M D ± > Mo° + Mw- In the rest frame of 
Z? ± we deal with 2-particle decay Z) ± — > DW^. In this 
frame the energy and momentum of , observed as qq 
dijet with effective mass Mw, are 



■M^-Ml 



A(M 2 D±1 M 2 



J w - 



2M 



2M 



D± 



(9) 



Denoting the W escape angle in D + rest frame relative 
to the direction of D + motion in the lab system by and 
c = cos 6, we have energy of W + in the lab system 

E w =i±(E w + cP ±P r ). 

Therefore, energies of dijets from Ws are distributed 
within the interval 

(£(-)=7±(^-/3±P r ), E (+) = 1± (E w + (3 ± p r )) . 

(10) 

Masses. The end point values Er±\ give two 
equations for evaluation of masses D ± and D° . 

(p r ) 2 and 



llM 2 



w 



In particular, 

Et_\ +Ei + \ — 2-f±E w . Therefore at large enough elec- 
tron energy (at 7 » 1), M 2 D± a E^M^/lE^E^}, 



Ml 



M w 



E P — E, 



(-) 



Ei 



(+) 



E e 



M D± . At finite 7 the 



exact equations are more complex. 

The accuracy of this procedure is determined by the 
accuracy of measurement of dijet energy together with 
its effective mass and by a width of D ± (if the latter is 
large). In particular, at so — the decay — > DW^- 
width is 



a 



r\3 



2 sm z W M 2 



(11) 



The T/M D ± ratio is below 0.1 at M D ± < 500 GeV. 

The distribution of these dijets in energy is uniform, 
dN(E) cx dE since there is no correlation between escape 
angle of W in the rest frame of D and production angle 
of Z? ± . When the width of D ± is not small, this dis- 
tribution become non-uniform near the end points. The 
measuring of fine structure of this distribution near the 
end point will give, at least roughly, the total D ± width. 

If M D ± >Mo°+Mw, the only decay channel is 
D + — > DW* + . All discussed above results are valid for 
each separate value of dijet effective mass M* , with evi- 
dent change in all equations Mw — > M* . The energy and 
M* distributions for each pair of dijets are independent 
from each other. The masses of Mo± and Mo are eval- 
uated in this case for each measured value of M* even 
with the best accuracy than in the previous case since in 
this case the proper width of D is low enough. 
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Spin. After evaluation of M D ±, the cross section of 
e+e- — > D + D~ process is calculated precisely for each 
s_d. The cross section for sd = j6|) is more than four 
times less than that for = 1/2 (JSJ) - This big difference 
allows to make a definite conclusion about spin of D- 
particles. 

Other properties. At sd = 1/2 the spins of D+ 
and D~ are correlated with longitudinal polarization of 
colliding electrons. In each event e + e~ — > D + D~ — > 
DDjj£ + v's we know the sign of dijet charge W = qq ■ 
It allows to study the charge and polarization asymme- 
tries for accessing of more detail properties of Z?-particles 
(e.g. ratio of D + D~ r y to D + D~Z couplings). 



III. THREE TYPES OF D-PARTICLES, D, D ± 
AND D A 

In the IDM together with D and D+ the one more neu- 
tral scalar particle D A exists (with mass M d a > Mo) 
0. Here CP parities of D or D A cannot be defined 
separately since they do not interact with fermions, 
but their relative parity is fixed, they have opposite 
CP-parities. Complete set of interactions in this case 
is given in (fT]). Similar particle can also exist in some 
models with sd — 1/2. So, we discuss now the case when 
additional neutral D-particle is D A . In this discussion 
we neglect case of approximate degeneracy between D 
and D A . 

• If M d a < M D i . 

1. 1. The lowest energy threshold for Z?-particle pro- 
duction has the process with decay D A — > DZ 

e+e~ -> DD A -> DDZ. (12) 

Instead of ((H), D A energy, 7-factor and velocity are 



4£ e 2 + Ml A - Ml 
A E e . 



1A 



E A 

I h 

E A r _A(4ElMl A ,Ml) (13) 

j PA 



M 



D A 



4E e E A 



In the IDM (at s£> = 0) the cross section is of the same 
order of value as a(e + e~ — > ^ + /i~): 

The signature of this process is similar to that given 
by eq. Q: 



One qq dijet or e + e~ or pair with iden- 

tical effective mass and energy distributions + 
nothing and with large missing E±. 



(15) 



effective mass equal to Mz- Energy distribution of this 
Z is given by equations similar to . (fTU|) . End points of 
this distribution allow to evaluate masses Mr> and M d a . 

If M d a — Md < Mz, the observable final state is Z* , 
which is seen as hadronic dijet or e + e~ or with 
identical spectra of effective mass. For each value of this 
effective mass M* , the energy distribution of this Z* is 
given by equations similar to ([51). (fTU|) . End points of this 
distribution allow to evaluate masses Mjj and M d a . 

2. After the study of process (fT2")) one must to study 
process (J^i) and cascade reactions 



D+D-^ 

DDW+W- , (a) 

D A W ± DW^^- DDW + W~ Z, (b) 

D A W+D A W- -> DDW+W- ZZ. (c) 



(16) 



If M d a — Mjj > Mz, the observable final state is Z, 
which is seen as hadronic dijet or e + e~ or with 



The decay D+ — > D A W+ is described by the same equa- 
tion as the decay D+- — > DW+. Its probability is lower 
than that for decay D+ — > DW+ due to smaller final 
phase space volume. Therefore, a(a) > a(b) > <r{c). 

The signature of the process (fTBTa)). just as the pro- 
cess iflBTb). (c)) for invisible decays of Z, is given by 
©. Each decay D+ -> DW+ and D ± -> D A W+ is de- 
scribed by identical equations, the only new point is that 
the end points E/_\ and for the energy distribution 
of Ws from decay D+ — > DW+ are given by equations 
(fTTJj) while the end points E^s and E A + ^ for energy dis- 
tribution of VT"'s from decay D+ —> D A W+ are given by 
the same equations but with the change Mo to M d a . It 
is evident that £7(_) < E A _^ < E A +) < E( + y Therefore 
in this case the same procedure as in sect. |TT] allows to 
obtain masses M D ± and Md (cross check for measuring 
of Mo)- 

To evaluate M d a from this reaction let us remind that 
for each type of decay the energy distribution of W in 
the lab system is uniform. The energy distribution of W 
in the lab system is the sum of two uniform distributions 
with the described above end points. The end points 
E A _y EA \ are marked by the steps in the density of event 
energy distribution. They can be used for new evaluation 
of M d a . 

The distribution of dijets in the effective mass can be 
different. If all masses are peaked around Mw, it means 
that M D ± - M D > M D ± - M d a > M w . If there are 
dijets with effective mass M\y and those with lower ef- 
fective mass, the former appear from D+ — > DW+ de- 
cay and the latter - from D+ — > D A W+ decay. If the 
effective masses of all dijets are below Mw, we have 
M w > M D ± - M D > M D ± - M d a . In this case the 
above mentioned steps in the dijet energy distribution at 
each M* will be added to steps in the distribution in M* . 

The signature of processes (TToT b)) and ([ToT c)) for vis- 
ible decays of Z is similar to with adding of dijet or 
£+£- pairs which represent Z or Z* . If one can distin- 
guish dijets from Z and those from W , the energy distri- 
bution of W in these processes can be used to enhance 
data massive for evaluation of masses, discussed above. 
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• M d a > M D ± , M d a + Mq < 2E e , the analysis of 
sect. [IT] is valid completely for the final states with signa- 
ture ((7]) - reaction e + e~ — > D + D~ . 

The second series of processes is ^jp) with two channels 
of D A decay and different signatures 

f +„^nD^l DDZ ' (a) (17) 

e e ^ruu -^<y DD ± W T -+DDW+W-. (b) (U > 

The process (fTTT a^ is the process (TL2")) . It can be ana- 
lyzed just as it was discussed earlier. The process (fI?Tb)) 
is a cascade process. It can be eliminated from mass anal- 
ysis of process ([3]) using the fact that in difference with 
the process ([3]) the observable decay products of this pro- 
cess move typically in one hemisphere. 

Note that in the case when M d a + Mp < 2M D ± 
the process ([T7|) has lower energy threshold than fl2Ji). 
The operations in suitable energy interval allows to find 
masses Mu and M d a but meet difficulties in evaluation 
of Md± . This problem can be solved by the increasing 
of beam energy for observation of process 

IV. SUMMARY 

We present the simple and robust method for discov- 
ery candidates for DM particles and evaluation of their 
masses and spins at ILC/CLIC. The same analysis can 
be applied to the case when the set of Z?-particles with 
RI < E e contains one additional neutral particle D\ with 
the same CP as D. 

These particles will be discovered via observation of 
processes with signature (|7J), ([15]) and with cross section 



of the order of a(e + e~ —> which is huge for LC. 

The masses of these particles will be obtained via measur- 
ing the details of the energy distribution of dijets (repre- 
senting W ± or Z) near the end points. The cross section 
measurements of processes with signature ([7]), (fl~5j) and 
similar signature for the derivative processes with cas- 
cade decay allow to determine the spin so of considered 
candidate for DM particle by comparison with simple SM 
calculation (the cross sections for so — 1/2 is approxi- 
mately 4 times larger than that for so = 0). 



One of proposed processes (J^i) was considered earlier 
in respect of discovery of neutralino as DMP, etc. (see 
e.g. [5|]). However, I never saw such approach for simul- 
taneous evaluation of masses and spins of D-particles 
irrespective to details of model. The advantages of 
presented approach arc following: 

1. The cross section of each suggested process is a 
substantial part of the total cross section of e + e~ 
annihilation at considered energy (typically up to tens 
percents). 

2. The signature is clear, the background is very small 
(typically, ~ 1 % from the observable effect). 

3. Simple kinematics allows to extract reliably quantities 
under interest from the data. 
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Serbo for discussions. The work was supported by grants 
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